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Color Glass Condensate and its relation to HERA physics 

Edmond lancu^ 

'^Institut de Physique Theorique de Saclay, F-91191 Gif-sur-Yvette, France 

I give a brief overview of the effective theory for the Color Glass Condensate, which is the high-density gluonic 
matter which controls high-energy scattering in QCD in the vicinity of the unitarity limit. 1 concentrate on 
fundamental phenomena, like gluon saturation, unitarization, and geometric scaling, and the way how these are 
encoded in the formalism. I emphasize the importance of the next-to-leading order corrections, especially the 
running of the coupling, for both conceptual and phenomenological issues. I survey the implications of the CGC 
theory for the IfERA physics and its phenomenological applications based on saturation models. 



1. Gluons at HERA 

The essential observation at the basis of the 
recent theoretical progress in the physics of 
hadronic interactions at high energy is the fact 
that high-energy QCD is the realm of high par- 
ton (gluon) densities and hence it can be stud- 
ied from first principles, via weak coupHng tech- 
niques. Anticipated by theoretical developments 
Hke the BFKL equation [1] and the GLR mech- 
anism [2|3j for gluon saturation, this observation 
has found its first major experimental foundation 
in the HERA data for electron-proton deep in- 
elastic scattering (DIS) at small-x. As visible in 
the leftmost figure in Fig. [TJ the gluon distri- 
bution xG{x, Q^) rises very fast when decreasing 
Bjorken-a; at fixed — roughly, as a power 
with A ~ 0.2 ~ 0.3. The physical interpretation 
of such results is most transparent in the pro- 
ton infinite momentum frame, where xG(x,Q^) 
is simply the number of the gluons in the proton 
wavefunction which are localized within an area 
Aa:_L ~ l/Q^ in the transverse plane and carry 
a fraction x ~ k^/Pz of the proton longitudinal 
momentum. 

Thus, without any theoretical prejudice, the 
HERA data suggest the physical picture illus- 
trated in the right hand side of Fig. [H which 
shows the distribution of partons in the transverse 
plane as a function of the kinematical variables for 
DIS in logarithmic units: InQ^ and Y = ln(l/a;). 
The number of partons increases both with in- 



creasing and with decreasing x, but whereas 
in the first case (increasing Q^) the transverse 
area ~ l/Q^ occupied by every parton decreases 
very fast and more than compensates for the in- 
crease in their number — so, the proton is driven 
towards a regime which is more and more dilute 
— , in the second case (decreasing x) the partons 
produced by the evolution have roughly the same 
transverse area, hence their density is necessarily 
increasing. 

Accordingly, the DGLAP equation [4] which 
describes the evolution with increasing is nat- 
urally linear, and also local in Q^. By contrast, 
the BFKL equation, which is the linear equation 
originally proposed pL, to describe the evolution 
with increasing energy, is non-local in transverse 
space and should be merely regarded as a linear 
approximation to more general evolution equa- 
tions which are non-linear, i.e., which account 
for the interactions among the partons within 
the wavefunction. The non-linear effects are ex- 
pected to become important in the region denoted 
as 'saturation' in Fig. [1] and in the approach to- 
wards it from the dilute region at large . 

Mainly because of its complexity, the high- 
energy evolution in QCD is not as precisely 
known as the corresponding evolution with . 
Still, the intense theoretical efforts over the last 
years led to important conceptual clarifications 
and to new, more powerful, formalisms — among 
which, the effective theory for the Color Glass 
Condensate (CGC) |5l6l7l8l9ll()j — , which en- 



1 



2 



E. lancu 




Figure 1. Left: The 1/x-evolution of the gluon, sea quark, and valence quark distributions for = 10 
GeV^, as measured at HERA (combined hts to the HI and ZEUS data). Note that the gluon and sea 
quark distributions have been reduced by a factor of 20 to fit inside the figure. Right: The 'phase- 
diagram' for QCD evolution as inferred from the HERA data; each colored blob represents a parton with 
transverse area Ax± 1/ and longitudinal momentum = xPz . 



compass the non-linear dynamics in high-energy 
QCD to lowest order in and allow for a unified 
picture of various high-energy phenomena rang- 
ing from DIS to heavy-ion, or proton-proton, col- 
lisions, and to cosmic rays. 

These developments provide a natural explana- 
tion for a variety of remarkable phenomena ob- 
served in the current experiments, like the 'ge- 
ometric scaling' in the HERA data at small x 
[11112] and the particle production at forward ra- 
pidities in deuteron-gold collisions at RHIC [T3| . 
Moreover, they have potentially interesting pre- 
dictions for the physics at LHC. It is our purpose 
in what follows to provide a brief, pedagogical, 
introduction to such new ideas, with emphasis on 
the physical picture and its consequences for deep 
inelastic scattering at high energy. 

2. DIS in the dipole frame 

At small x, DIS is most conveniently com- 
puted by using the dipole factorization (see, e.g., 
Refs. [TU] for more details and references). The 



small-x quark to which couple the virtual photon 
is typically a 'sea' quark produced at the very end 
of a gluon cascade. It is then convenient to dis- 
entangle the electromagnetic process "f*q, which 
involves this 'last' emitted quark, from the QCD 
evolution in the proton, which involves mostly 
gluons. This can be done via a Lorentz boost 
to the 'dipole frame' in which the struck quark 
appears as an excitation of the virtual photon, 
rather than of the proton. In this frame, the pro- 
ton still carries most of the total energy, while the 
virtual photon has just enough energy to dissoci- 
ate long before the scattering into a 'color dipole' 
(a qq pair in a color singlet state), which then 
scatters off the gluon fields in the proton. This 
leads to the following factorization: 

ar'p{x,Q^) ^ J^dzJ dV |*^(z,r;Q2)|2 

X crdipoie(a;,r) (1) 

where |^'-y(z, r; Q^)p is the probability for the 
7* — > qq dissociation (r is the dipole trans- 
verse size and z the longitudinal fraction of the 
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quark), and (Tdipoiol^;, is the total cross-section 
for dipole-proton scattering and represents the 
hadronic part of DIS. At high energy, the latter 
can be computed in the eikonal approximation as 



CTdipoic(a;,r) = 2 (Pb T{x,r,b) 



(2) 



where T(x, r, b) is the forward scattering ampli- 
tude for a dipole with size r and impact parame- 
ter b. This is the quantity that we shall focus on. 
The unitarity of the 5-matrix requires T < 1, 
with the upper limit T = I corresponding to to- 
tal absorbtion, or 'black disk limit'. 

But the unitarity constraint can be easily vio- 
lated by an incomplete calculation, as we demon- 
strate now on the example of lowest-order (LO) 
perturbation theory. To that order, T(r, b, Y) in- 
volves the exchange of two gluons between the 
dipole and the target. Each exchanged gluon 
brings a contribution gf^r ■ Ea, where Ea is the 
color electric field in the target. Thus, T ^ 
g'^r'^{Ea ■ Ea)x, where the expectation value in 
the r.h.s. is recognized as the number of gluons 
per unit transverse area, as measured by a probe 
with transverse resolution ^ 1 /r^ : 

2 xG{x,l/r'^) 



T{x,r, b) 



7ri?2 



as n{x, Q ) . (3) 



Above, we identified the gluon occupation num- 
ber: n{x,Q^) = [number of gluons xG{x,Q'^)] 
times [the area occupied by each gluon] di- 

vided by [the proton transverse area ttR'^]. 

Eq. ^ applies so long as T <C 1 and shows that 
weak scattering (or 'color transparency') corre- 
sponds to low gluon occupancy n <^ 1/q?s. But if 
naively extrapolated to very small values of x, this 
formula leads to unitarity violations : T would 
eventually become larger than one ! Before this 
happens, however, new physical phenomena are 
expected to come into play and restore unitarity. 
As we shall see, these are non-linear phenomena, 
and are of two types: (i) multiple scattering, i.e., 
the exchange of more than two gluons between 
the dipole and the target, and (ii) gluon satu- 
ration, i.e., non-linear effects in the proton wave- 
function which tame the rise of the gluon distri- 
bution at small x. 

Eq. ^ also provides a criterion for the onset 
of unitarity corrections: these become important 



when T{x,r) ~ 1 or n{x,Q^) ~ l/a^. This con- 
dition can be understood as follows: to have non- 
linear phenomena, the gluons in the wavefunction 
must be numerous enough (which requires small 
x) and large enough (meaning low transverse mo- 
menta Q^) in order to strongly overlap with each 
other, by a factor n ^ l/ofg ^ 1 which is large 
enough to compensate for the smallness of the 
coupling. When this happens, the gluon mutual 
interactions ~ nag become of 0(1). 

The condition n{x,Q^) ~ 1/as can be solved 
for the saturation momentum, which is the value 
of the transverse momentum below which satura- 
tion effects are expected to be important in the 
gluon distribution. One thus finds 



1 



(4) 



which grows with the energy as a power of 1/x, 
since so does the gluon distribution before reach- 
ing saturation. In logarithmic units, with Y = 
ln(l/a;), the saturation line liiQ^(Y) — XY is 
therefore a straight line, as illustrated in the right 
hand side of Fig. [TJ This is the borderhne be- 
tween the dilute regime at high transverse mo- 
menta k± Qs{Y), where one expects the stan- 
dard perturbation theory to apply, and a high- 
density region at low momenta k± ^ Qs(Y), 
where physics is non-linear. In fact, as we shall 
argue below, at high energy the effects of satura- 
tion can extend up to very high values of k± , well 
above the saturation line. 

3. BFKL evolution 

Within perturbative QCD, the emission of 
small-a; gluons is amplified by the infrared sen- 
sitivity of the bremsstrahlung process, whose it- 
eration leads to the BFKL evolution (at least, for 
not too high energies). Fig. [2] shows the emission 
of a gluon which carries a fraction x = k^/pz of 
the longitudinal momentum of its parent quark. 
When a; <C 1, the differential probability for this 
emission can be estimated as 



dPR 



UsCf d^k± dx 
27r2 kl ~ 



(5) 



which is singular as a; ^ 0. Introducing the rapid- 



4 



E. lancu 




0000000000-' = xp^ 



Figure 2. Bremsstrahlung at lowest order. 

ity Y — 111(1/2:), and hence dY — dx/x, Eq. ^ 
shows that there is a probability of 0{as) to emit 
one gluon per unit rapidity. The same would hold 
for the emission of a soft photon from an electron 
in QED. However, unlike the photon, the child 
gluon is itself charged with 'colour', so it can fur- 
ther emit an even softer gluon, with longitudinal 
fraction xi = Qz/kz 1. When the rapidity is 
large, agY ^ 1, such successive emissions lead to 
the formation of gluon cascades, in which the glu- 
ons are ordered in rapidity and which dominate 
the small-x part of the hadron wavefunction. 

So long as the density is not too high, the glu- 
ons do not interact with each other and the evo- 
lution remains linear : when further increasing 
the rapidity in one more step (Y ^ Y + dY), 
the gluons created in the previous steps incoher- 
ently act as color sources for the emission of a 
new gluon. This picture leads to the following, 
schematic, evolution equation 



dn 
dY 



n(Y) cx e' 



(6) 



which predicts the exponential rise of n with Y . 
There is an additional feature of the high en- 
ergy evolution which needs to be emphasized: 
the gluon emission vertex is non-local in trans- 
verse momentum (the transverse momentum of 
the daughter gluon is generally different from that 
of its parent) , but this non-locality is quite weak 
and can be described as diffusion in the logarith- 
mic momentum variable p = InQ^. That is, a 
better version of Eq. ^ reads 



dn{Y,p) 
dY 



(7) 



where x is the diffusion coefficient. This is an 
oversimplified version of the BFKL (Balitsky- 
Fadin-Kuracv-Lipatov) equation [T] which cap- 
tures the main features of this evolution: the 



unstable growth of the gluon distribution and 
the diffusion in transverse momentum. Both fea- 
tures leads to difficulties in the high-energy limit: 

(i) the unlimited growth of n entails unitarity 
violations, as discussed in the previous section; 

(ii) the BFKL evolution explores the transverse 
phase-space in a diffusive way, so its solution 
n(y, p) receives significant contributions from all 
the points p' such that \p' — /o| < y/xctsY. Hence, 
even if p is hard, the momenta p' contributing 
to the solution can be considerably softer; with 
increasing Y, there is a larger and larger part of 
the total result which is generated from soft mo- 
menta (on the left of the saturation line), where 
linear evolution, or even perturbation theory, are 
bound to fail. One knows by now that both fea- 
tures — the exponential growth of n with Y and 
the BFKL diffusion — are considerably tempered 
by NLO effects |14|15j . like the running of the 
coupling or the constraint of energy conservation. 
But the basic fact that the gluon density increases 
very fast with Y is expected to remain true (to 
all orders in as) so long as one neglects the non- 
linear effects in the evolution. 

4. JIMWLK evolution and the CGC 

Non-linear effects appear because gluons carry 
colour charge, so they can interact with each 
other (even when separated in rapidity) by ex- 
changing gluons in the i-channel, as illustrated 
in Fig. [3) These interactions are amplified by the 
gluon density and thus they should become more 
and more important when increasing the energy. 

Back in 1983, L. Gribov, Levin and Ryskin [2] 
(see also Ref. [3, ) suggested that gluon saturation 
should proceed via 2^1 'gluon recombination', 
which is a process of order a^n^ (cf. Fig. [3]). One 
can heuristically take this into account by adding 
a non-linear term to the r.h.s. of Eq. ^ : 

— ~ asd1n + asU - aln^ . (8) 

Clearly, this non-linear equation has a fixed point 
JT-sat ~ l/tts at high energy. That is, when n is as 
high as 1/as, the emission processes (responsible 
for the BFKL growth) are precisely compensated 
by the recombination ones, and then the gluon 
occupation factor saturates at a fixed value. 
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Figure 3. DIS in the presence of BFKL evolution, saturation and multiple scattering. 



Twenty years later, we know that the ac- 
tual mechanism for gluon saturation in QCD is 
more subtle than just gluon recombination and 
that its mathematical description is consider- 
ably more involved than suggested by Eq. 
This mechanism, as encoded in the effective 
theory for the CGC and its central equation, 
the JIMWLK equation (Jalilian-Marian, lancu, 
McLerran, Weigert, Leonidov, and Kovner) [617] . 
is the saturation of the gluon emission rate due 
to high density effects : At high density, the glu- 
ons are not independent color sources, rather they 
are correlated with each other in such a way to 
ensure color neutrality [16117118] over a distance 
Ax^ ^ 1/Qs- Accordingly, the soft gluons with 
k± ^ Qs Sire coherently emitted from a quasi- 
neutral gluon distribution, and then the emission 
rate dn/dY saturates at a constant value of 0(1). 
Thus, in the regime that we call 'saturation', the 
gluon occupation factor keeps growing, but only 
linearly in Y (i.e., as a logarithm of the energy) 
|19ll6j . Schematically: 



^ = Xin) « (9) 
[l ifn>l/a, 

where x("-) is a non-linear function with the lim- 
iting behaviours displayed above. The transition 
between the two regimes is smooth and it occurs 
around the saturation line, i.e., at transverse mo- 



menta k± ^ Qs{Y), where Qs(Y) is an increasing 
function of Y which is determined by the the- 
ory. As we shall later explain, this rise is roughly 
consistent with the power-like increase with 1/x 
predicted by Eq. (jH). 

Eq. ^ is not yet the JIMWLK equation, but 
only a mean field approximation to it: in reality, 
one cannot write down a closed equation for the 
2-point function n{Y) = {Ea ■ Ea)Y, but only 
an infinite hierarchy for the A^-point correlations 
{A{l)A{2) ■ ■ ■ A{N))y of the gluon fields. This is 
so since the A^-point functions couple under the 
evolution via the non-linear effects. In the CGC 
formalism, these correlations are encoded into the 
weight function Wy [A] — a functional probability 
density for the color field configurations: 

{A{1) A{2) ■ ■ ■ A{N))y = 

= J D[A] Wy[A] A{1) A{2) ■ ■ ■ A{N) . (10) 



The average in Eq. pO]) is similar to the 'aver- 
age over disorder' that is usually performed in 
the study of amorphous materials, like glasses: 
the various target configurations scatter indepen- 
dently with the incoming projectile (indeed, their 
internal dynamics is 'frozen' over the character- 
istic time scale for scattering, by Lorentz time 
dilation), and the physical scattering amplitude 
is finally obtained by summing the contributions 
from all such configurations, with weight function 
M/V[^]- This explains the concept of 'glass' in 
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the 'Color Glass Condensate'. The 'color' refers, 
of course, to the gluon color charge. Finally, the 
'condensate' stays for the coherent state made by 
the gluons at saturation: this state has a large 
occupation number n ~ l/tts ^ 1, as typical for 
a Bose condensate. 

The JIMWLK equation [6|7j is a functional 
differential equation describing the evolution of 
Wy[A] with Y. Via Eq. this functional 

equation generates an infinite hierarchy of ordi- 
nary evolution equations for the TV-point func- 
tions of A, as anticipated. To describe physical 
observables, these correlations must be gauge- 
invariant. At high energy, a particularly conve- 
nient set of gauge-invariant correlation functions 
is obtained by taking products of Wilson lines 
traced over the color indices. 

A 'Wilson line' describes the scattering be- 
tween a high energy parton and a gauge back- 
ground field in the eikonal approximation: the 
parton preserves a straight-line trajectory while 
moving through the field. Hence, the product of 
N Wilson lines describes the S'-matrix for N par- 
tons propagating in the background field A. After 
also averaging over A, as in Eq. IjlOp. one finally 
obtains the S'-matrix for the eikonal scattering 
between the partonic system and the hadron (the 
'CGC'). For instance, the average S-matrix for 
dipole scattering, as relevant to DIS (cf. Fig. [3]), 
is computed as 

{S{r,b))Y = J D[A] Wy[A] ^tr(yjy^) (11) 

where is the Wilson line for the quark with 
transverse coordinate x, i.e., 

V^x) ^ P exp (^ig J dx-A+{x~,x)t''^ ,(12) 

Vy is the corresponding operator for the antiquark 
at y, and the dipole size and impact parame- 
ter are given by r = a; — y, b ~ {x + y)/2. The 
(average) forward scattering amplitude, which is 
the quantity which determines the dipole cross- 
section ([21), is then obtained as (T) = I — (S). 

In this description, the unitarity corrections 
are explicit: the multiple scattering is encoded 
in the Wilson lines and the gluon saturation in 



the weight functional PVV[^]- Hence, no surpris- 
ingly, the corresponding evolution equations, as 
deduced from the JIMWLK equation, have the 
property to manifestly preserve unitarity. These 
equations form an infinite hierarchy which was 
originally derived (within a different formalism) 
by Balitsky [8^. The first equation in this hierar- 
chy reads (with transverse coordinates omitted) 

9y(r) = a,{T) - as{T^). (13) 

This is not a closed equation: the amplitude 
(T) for the scattering of one dipole is related 
to the amplitude (T^) for two dipoles. A closed 
equation, known as the Balitsky-Kovchegov (BK) 
equation [9j, can be obtained in a mean field ap- 
proximation which assumes factorization: (T^) » 
(T)(r). After restoring the transverse coordi- 
nates in the diffusion approximation, the BK 
reads (with p = \n{l/r'^) and T = (T)) 

dyTiY, p) ^ a^dlT + a,T - a,T\ (14) 

As anticipated, unitarity (T < 1) is manifest on 
this equation which has T = 1 as a fixed point 
at high energy. Formally, the BK equation repre- 
sents the large- TVc limit of the Balitsky- JIMWLK 
hierarchy. Recent numerical studies demonstrate 
that this mean field aproximation works better 
than expected: for Nc = 3, the differences be- 
tween the BK and JIMWLK predictions for the 
dipole amplitude are less than 1% [20^. These 
properties, together with the relative simplicity 
of Eq. p4p , make this equation a very convenient 
tool for studies of saturation and unitarity. Some 
of its physical consequences will be described in 
the next sections. 

Let us conclude this section on the general for- 
malism with a few additional remarks: 

(i) The original motivation for the CGC the- 
ory [5] came from the physics of ultrarelativistic 
heavy ion collisions, where this theory is meant 
to describe the initial states of the incoming 
nuclei, prior to collision. In that context, the 
gluon density is large already at moderate en- 
ergies, due to the existence of many 'tree-level' 
color sources: the 3A valence quarks, with A the 
atomic number. The CGC theory produced some 
interesting predictions for particle production in 
nucleus-nucleus collisions, which have been since 
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confirmed by the experimental results at RHIC 
(see the review articles 10,21 22. 23^ and Refs. 
therein). With the advent of LHC, the CGC 
theory should find a vaste field of applications 
to both nucleus-nucleus and proton-proton colli- 
sions [24] , 

(ii) The BK equation in the diffusion approx- 
imation, as written down in Eq. p4p . turns out 
to be the same as the FKPP equation which de- 
scribes the mean field limit (corresponding to very 
large occupation numbers) of a classical stochas- 
tic process known as reaction- diffusion. This 
process can be briefly described as follows [40] : 
'molecules' of type A which are located at the 
sites of an infinite, one-dimensional, lattice can 
locally split {A AA) or merge with each other 
{AA A); also, a molecule can diffuse to the 
adjacent sites. The correspondence between BK 
and FKPP, originally noticed in Ref. [39], sheds 
new light on the physics of geometric scaling and 
helps clarifying the limitations of the mean field 
approximation (see below). 

(iii) The central equations of the CGC for- 
malism, so like the Balitsky-JIMWLK hierarchy, 
or the factorization formula ([TT]) for DIS together 
with the corresponding ones for proton-nucleus 
[25126] . or nucleus-nucleus [5127] , collisions, are 
known so far only to leading order (LO) in per- 
turbation theory — an approximation which by 
itself involves an infinite resummation of pertur- 
bative contributions where the powers of g are 
accompanied by appropriate powers of y = Ins, 
or of the strong gauge fields A ^ 1/g, so hke in 
Eq. ITS]) . However, the higher-order effects, and 
especially the running of the coupling, turn out to 
be extremely important — whenever known, they 
dramatically affect the predictions of the LO the- 
ory, as we shall shortly see. Fortunately, there is 
an ongoing effort towards the inclusion of higher- 
order corrections, which so far has led to an im- 
proved version of the BK equation containing run- 
ning couphng effects [28129] . 

(iv) Even at leading order, the equations pre- 
viously mentioned (Balitsky-JIMWLK and BK) 
are still incomplete: they neglect the effects of 
gluon number fluctuations (or 'Pomeron loops'), 
i.e., the correlations associated with the fact that 
some of the gluons produced by the high-energy 



evolution have a common ancestor [33|34|35] . 

However, in practice, this is not a serious draw- 
back, since these correlations are anyway sup- 
pressed by the running of the coupling [36] . This 
will be further discussed in the next section. 

5. Saturation and Geometric scaling 

Due to its simplicity, the BK equation p4|) is 
well suited for both numerical and analytic stud- 
ies of the evolution towards unitarity and satu- 
ration. The corresponding solution TiY^p) is a 
front which with increasing Y propagates towards 
larger vales of p, in such a way that the position 
of the front coincides with the saturation momen- 
tum p,(r) EE \nQl{Y) [37138139] (see Fig. [g . Be- 
hind the front, i.e., for p < Ps(i^), the dipole scat- 
tering amplitude has reached the black disk limit 
T — \ (which, we recall, is a fixed point of the 
BK equation) and hence it cannot grow anymore. 
Ahead of the front {p ^ ps{Y)), the amplitude 
is still weak, T <^ 1, so the non-linear term in 
Eq. (|14p is unimportant and the amplitude can 
grow according to the linear, BFKL, evolution. 



1/2 




P(Y,) P(Y2) P 



Figure 4. The saturation front generated by the 
BK equation for two values of the rapidity. 

This argument suggests that the progression 
of the front with increasing Y is entirely driven 
by the linearized version of the BK equation — 
a conclusion which is confirmed by more rigor- 
ous mathematical arguments '33^. Hence, by 
solving the linear, BFKL, equation supplemented 
with a saturation boundary condition (namely. 
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T{Y,p) ~ 1 for p = ps{Y)), one can compute 
the positfon Ps{Y) of the front, which is the 
same as the saturation line, and also the shape 
of the amphtude ahead of the front. One thus 
finds [37_38 that the front progresses at constant 
speed, i.e., the saturation hne is a straight hne 
(befow, as = asNc/n) : 

p,{Y) XsY with As«4.88as. (15) 

Moreover, within a relatively large window at 
p > ps(Y), whose width is growing with Y, the 
amplitude depends upon p and Y only via the 
difference t = p- ps{Y) = /r^Q^,{Y)) : 

r(r, p) cx T e-'^=^ for 1 < r < , (16) 

with 7s ~ 0.63. This property is known as geo- 
metric scaling. This scaling also holds (trivially !) 
behind the front, since T is constant there: T = 1 
for p ^ Ps{y)- Geometrically, this means that, 
with increasing Y , the saturation front gets sim- 
ply translated towards larger values of p, but its 
shape remains unchanged: the front propagates 
like a traveling wave |39j . 

We conclude that the dipole amplitude shows 
geometric scaling for all values of p up to a maxi- 
mal value Pgoom 

(Y) ~ ps{Y) + Vx^IT, which for 
large Y can be significantly larger than Ps{Y). 
(The difference Pgeom — Ps Vy grows with Y 
via BFKL diffusion.) This means that the effects 
of saturation make themselves felt even at rela- 
tively large momenta 3> Ql{Y), where the 
scattering is weak, T ^ 1, and the gluon den- 
sity in the target is quite low. This considerably 
extends the phase-space where saturation is ex- 
pected to be important in the experiments. 

Via the factorization formula the scaling 
of the amplitude as a function of rQs(Y) im- 
plies a similar scaling for the DIS cross-section: 
<T.y*piY,Q^) « (t(t), with T = Q^/Q^^{Y). Such 
a scaling has been indeed identified in the HERA 
data, by Stasto, Golec-Biernat and Kwiecihski 
[TT] (see Fig. [5]). More recently, with the advent 
of more precise data for DIS diffraction at HERA, 
geometric scaling has been noticed in these data 
too [12]. The data also show violations of geo- 
metric scaling, which can be understood as con- 
sequences of the BFKL diffusion 42J and of the 
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Figure 5. Geometric scaling in the HERA data 
for (7-y*p at X < 0.01; t is the scaling variable, 

T ^ QVQiiY) m- 

quark masses 03] . These phenomenological stud- 
ies will be further described in the next section. 

But the results at HERA also show that the 
rise of the gluon distribution with 1/x is much 
slower than predicted by the LO BFKL analy- 
sis: when described in terms of saturation, they 
require a saturation exponent As — 0.2 0.3 
whereas Eq. yields As — 1 for as — 0.2. 
This discrepancy is solved by the NLO calcula- 
tion of the saturation exponent [41| which pre- 
dicts indeed As — 0.3. This large difference 
between the LO and the NLO results for As is 
largely explained by the running of the coupling, 
whose consequences are crucial for the physics of 
saturation. The relevant value of the coupling 
is that corresponding to the saturation momen- 
tum, as{Qs[Y)) cx l/ln[Q2(y)/A|(.j3], which de- 
creases with Y . Hence, for sufficiently high en- 
ergy, all the other NLO corrections (like the NLO 
effects in the BFKL kernel jl4ll5) ) are suppressed 
by ois{Qs{Y)) <C 1, so that the correct theory for 
saturation (at least for asymptotically large Y) is 
the LO theory extended to running coupling. 

The main effect of the running of the coupling 



Color Glass Condensate and HERA 



9 



is to considerably slow down the evolution, with 
some dramatic consequences: (i) rather than 
growing exponentially with the saturation mo- 
mentum grows only like exp \/Y [37138] : (ii) the 
window pgoom — Ps for geometric scaling outside 
saturation grows very slowly, as Y^^^ [41]; (iii) 
for a nuclear target, the dependence of Qs upon 
the atomic number A is strongly reduced by the 
running of coupling — for sufficiently high energy, 
there should be no difference between the satura- 
tion scale of a nucleus and that of the proton [30] . 
This last feature has intriguing consequences for 
particle production in proton-nucleus collisions 
at very high energies, perhaps at LHC [31^ . 

Another important, and rather surprising, con- 
sequence of the running of the coupling [36] is to 
improve the applicability of the mean field ap- 
proximation, that is, the BK equation. We have 
already mentioned, at the end of Sect. 4, that 
the basic equations of the CGC theory (Balitsky- 
JIMWLK and BK) ignore the correlations in- 
duced via gluon number fluctuations. These cor- 
relations refer to the fact that the gluons pro- 
duced by the evolution can have common ances- 
tors. Or it turns out that, for a fixed coupling, 
these correlations do significantly affect the pic- 
ture of saturation [32I33I34I35] : Although they 
are mostly produced in the dilute regime, i.e., 
in the tail of the gluon distribution at high Q^, 
these correlations are rapidly amplified by the 
BFKL evolution, so they eventually influence the 
approach towards saturation. (In terms of dia- 
grams, the evolution with both fluctuations and 
saturation contains Pomeron loops.) Accordingly, 
this evolution becomes stochastic and the satura- 
tion scale itself becomes a random variable, whose 
dispersion increases with Y : with a fixed cou- 
pling, this rise is so fast that already for mod- 
erate values of Y it completely washes out the 
mean field picture (and, in particular, the prop- 
erty of geometric scaling). Direct calculations 
of such effects in QCD are extremely difhcult 
(the complete theory for QCD evolution with 
Pomeron loops is still lacking; see, however, Refs. 
[35l44l45l4HI4Tt48] V but the effects of the fluc- 
tuations can be appreciated from the experience 
with the reaction-diffusion problem in statistical 
physics, and also with some QCD-inspired mod- 



els which allow for explicit numerical calculations 
and belong to the universality class of reaction- 
diffusion |49|50|51|52l53] . But these studies also 
allow for the inclusion of a running coupling, and 
the effects of that turn out to be dramatic [36] 
: the fluctuations are strongly suppressed up to 
the highest values of Y of interest, so that the 
complete, stochastic, evolution gives essentially 
the same results as its mean field approximation. 
This is to be attributed to the fact that, with run- 
ning coupling, the saturation front has a different 
shape (due to the shrinking of the window for 
geometric scaling), which disfavors fluctuations 
[36|. We thus conclude that the actual evolution 
in QCD at high energy is not in the universality 
class of reaction-diffusion, and this is somehow 
fortunate as it allows us to rely on mean field ap- 
proximations like the BK equation (with running 
coupling, of course). 

6. Saturation models 

Although fully consistent with the gross fea- 
tures of the experimental results at both HERA 
and RHIC, the current formalism for high-energy 
evolution with saturation is not accurate enough 
to allow for a precise, parameter-free, descrip- 
tion of the small-x HERA data, which are known 
with high precision. The theoretical limitations 
refer to both perturbative and non-perturbative 
aspects: on the one hand, the NLO corrections 
(expected to be large) have not been system- 
atically implemented; on the other hand, the 
impact-parameter dependence of the scattering 
amplitudes goes beyond perturbation theory (es- 
pecially in the dilute region towards the periphery 
of the hadron disk). To cope with that, various 
"saturation models" have been proposed for the 
dipole cross-section, Eq. ([2]), which were inspired 
by theoretical ideas about saturation, or by re- 
sults from the BK equation, but which are also in- 
volving several (typically, 3 or 4) free parameters. 
Such models provided remarkably good descrip- 
tions of the relevant data at HERA and RHIC — 
even surprisingly good, given the simplicity of the 
models and the reduced number of free parame- 
ters. Some general remarks about these models: 
(i) The free parameters are fixed from fits to 
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Figure 6. Predictions of the 'CGC model' with 4 quarks (3 hght + 1 heavy). Left: the charm structure 
function. Right: the bottom structure function and the longitudinal structure function. From Rcf. 



the F2 data alone; all the other results emerge as 
predictions, and they provide a reasonably good 
description of the ensemble of the HERA data at 
X < 0.01, including the longitudinal (Fl), diffrac- 
tive (F^), and charm (^"2) structure functions, 
the virtual photon production of vector mesons 
J/^), and the deeply virtual Compton scat- 
tering (DVCS). (ii) The saturation models pro- 
vide natural explanations for important, qualita- 
tive, features of the data like geometric scaling, 
the turn-over in F2 at low x and low Q^, and 
the nearly constant diffractive-to-inclusive ratio 
(Tdiff/ctot at HERA, or the total multiplicity and 
the high-p_L suppression of particle production in 
forward d-Au collisions ('i?pA ratio') at RHIC. So 
far, there are no other compelling explanations 
which apply to the ensemble of these data. 

The first "saturation model", due to Golec- 
Biernat and Wiisthoff (GBW) played an 

important role towards the shift of paradigm in 
favour of saturation at HERA. The main virtue 
of that fit was in its simplicity: with a very simple 
functional form. 



Ql{x) = {x^lxfGeY\, 



(17) 



which interpolates between color transparency 
(cTdipoic c>c r^) for small dipole sizes and saturation 
(^dipoic — 27ri?^) for larger dipoles, and the tran- 
sition occurring at a 'critical' scale Tg ^ \/Qs(x) 
which decreases with 1/x (the real hallmark of 
saturation) , this model offered a good description 
of the early HERA data (at x < 0.01 and any Q^) 
with only 3 free parameters: the proton radius i?, 
the value xq where Qs — 1 GeV, and the satura- 
tion exponent A (the data favoured xq « lO"'^ 
and A « 0.3). Note that this model has exact 
geometric scaling built in, and in fact it was his 
success which inspired the search for this scaling 
in the data [TT]. 

However, the limitations of this model become 
obvious with the advent of more precise HERA 
data, and new, more sophisticated, models were 
then proposed to account for these data. The 
main improvements referred to a better inclusion 
of the effects of the perturbative QCD evolution 
(which in particular brought in violations of ge- 
ometric scaling), sometimes accompanied by a 
more complex treatment of the parameter-impact 
dependence. Some approaches [55|56|57j focused 
on improving the high-Q^ behaviour of the fit 
(i.e., the small-r behaviour of the dipole cross- 
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Figure 7. Comparison of the various fits in Ref. j6'0f to a subset of DIS data. Left: No saturation fits. 
FS2004 Regge dipole fit (dashed line) and (solid line) a fit of the same model to data in the restricted 
range 5 x 10^^ < a; < 10~^, extrapolated over the whole x-range x < 0.01. Right: Saturation fits. FS2004 
saturation fit (solid line) and the CGC dipole model (dot-dashed line) 



section), by adding in the DGLAP evolution. 
Some others have rather emphasized the BFKL 
physics and the transition from hnear to non- 
hnear dynamics [42l43j . 

In particular, the 'CGC model' in Ref. 
which is based on approximate solutions to the 
BK equation, has shown that both the BFKL 
value for the 'anomalous dimension' (the slope 
7 « 0.63 in Eq. ^6]) ) and the pattern of geomet- 
ric scaling violations predicted by the BFKL dif- 
fusion are consistent with the HERA data. This 
fit, which involves 3 light quarks and the same 3 
free parameters as the GBW fit, has also shown 
that the data prefer a smaller value for the sat- 
uration exponent, namely A « 0.25. This value 
has been further reduced, to A ~ 0.22, after also 
including the heavy charm quark in the fit [13]; 
this last analysis requires a somewhat larger value 
7 « 0.76 for the 'anomalous dimension'. Some 
predictions of [43] are summarized, together with 
the respective data at HERA, in Fig. [6l The 
CGC model has been recently extended to in- 
clude impact-parameter dependence |57|58j . but 
the respective fit favors an unusually small value 
7 « 0.46 for the 'anomalous dimension', which 
looks inconsistent with the BK dynamics. 

An alternative saturation model, 'FS2004', has 



been proposed in Refs. [59160161] . which is partic- 
ularly simple (and thus closer in spirit to the orig- 
inal GBW model), but also more fiexible, in the 
sense of including more free parameters. Interest- 
ingly, this fit has two versions (with and without 
saturation), and the version including saturation 
appears to be clearly favored |60j by the HERA 
data for F2, so long as all the data (including 
those at low Q^) are included in the fit. This is 
illustrated in Fig. [7] [60 , which also show that the 
quality of the FS2004 fit with saturation is similar 
to that of the CGC fit in Ref. 

One should also mention here the saturation 
models used in the context of RHIC, where 
the dipole cross-section (or, more precisely, its 
Fourier transform) enters the rate for forward 
particle production in deuteron-gold collisions. 
In that case, it is preferable to formulate the 
models directly in momentum space, to avoid 
numerical artifacts associated with the Fourier 
transform. The saturation models formulated in 
Refs. |62|i63j provide a reasonable description of 
all the relevant RHIC data; moreover, as shown in 
Ref. [M] , a slight modification of these models can 
also account for the HERA data, as expected on 
the basis of the universality of saturation physics. 

Returning to HERA physics, we shall conclude 
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Figure 8. Left: Schematic view of DIS diffraction. Right : Quasi-elastic diffraction ('Vgap — Y). 



this review with a brief discussion of DIS diffrac- 
tion, which is a particularly convenient labora- 
tory to test saturation. Indeed, the diffractive 
cross-section is controlled by large dipole sizes 
and hence it is particularly sensitive to our the- 
oretical ideas about unitarization. The theory 
is simpler for the case of a large rapidity gap 
Ygap ~ y, or small diffractive mass M"^ ~ Q^, 
in which case diffraction amounts to the elastic 
scattering of the qq dipole produced by the disso- 
ciation of 7* (see Fig. [8]). The respective cross- 
section is then evaluated as (compare to Eq. ([T])) 



game and cuts off the integral in Eq. ([T 



d26 



= Jdzd^r \^^{z,r;Q)\' (r(r, F))' . (18) 



The photon wavcfunction favors relatively small 
dipoles with r ^ 1/Q : 



i/Q^ 



But for small r, T{r,Y) cx r^, hence the inte- 
gral will be dominated by the size rg where the 
amplitude has a turn-over, due to unitarity cor- 
rections. If unitarization is to be associated with 
the soft, non-perturbative, physics (the prevail- 
ing viewpoint before the advent of saturation; 
see, e.g., [53]), then ^ 1/Aqcd, and diffrac- 
tion would be non-perturbative even when is 
hard ! However, for sufficiently large Y (small x), 
the (semi)hard saturation scale Qs{Y) enters the 



dgdiff 



(20) 



i/Q" 



In this scenario, (Xdiff oc Q'^{x) cx x~'^ scales like 
the 'hard Pomeron', that is, in the same way as 
the inclusive cross-section ([T]); hence the ratio 
fdiff / Ctot is approximately constant as a function 
of the energy. This prediction is to be contrasted 
to that of non-perturbative unitarization [65], 
where one rather expects the diffractive cross- 
section to rise twice as fast as the gluon distri- 
bution (since elastic scattering requires the ex- 
change of at least two giuons), which would im- 
ply cTdiff/crtot ^ x"'^ . It turns out that the HERA 
data favor the saturation scenario: the measured 
ratio tJdiff / CTtot is very flat as a function of the in- 
variant energy oc 1/x, as illustrated in Fig. [HI 
This figure also shows that this flatness is well re- 
produced by the saturation model of Ref. [55] . 
Notice that the integrand in Eq. (|20p rises as 
a double (hard) Pomeron, Ql{x) oc x~'^^, but 
one power of Qs(a;) is eventually compensated 
by the energy-dependence of the upper cutoff 
r1 — \/Ql{x), which is the critical size for the 
onset of unitarity in the framework of saturation. 
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